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Classical lattice simulations of the Abelian Higgs model are used to investigate small scale struc- 
ture and loop distributions in cosmic string networks. Use of the field theory ensures that the 
small-scale physics is captured correctly. The results confirm analytic predictions of Polchinski & 
Rocha Pi for the two-point correlation function of the string tangent vector, with a power law from 
length scales of order the string core width up to horizon scale with evidence to suggest that the 
small scale structure builds up from small scales. An analysis of the size distribution of string loops 
gives a very low number density, of order 1 per horizon volume, in contrast with Nambu-Goto sim- 
ulations. Further, our loop distribution function does not support the detailed analytic predictions 
for loop production derived by Dubath et al. 0). Better agreement to our data is found with a 
model based on loop fragmentation coupled with a constant rate of energy loss into massive 
radiation. Our results show a strong energy loss mechanism which allows the string network to scale 
without gravitational radiation, but which is not due to the production of string width loops. From 
evidence of small scale structure we argue a partial explanation for the scale separation problem of 
how energy in the very low frequency modes of the string network is transformed into the very high 
frequency modes of gauge and Higgs radiation. We propose a picture of string network evolution 
which reconciles the apparent differences between Nambu-Goto and field theory simulations. 

PACS numbers: 



I. INTRODUCTION 

Cosmic strings are line-like objects formed in the early 
universe (for reviews see They exist as solitons 

in theories with spontaneously broken symmetries if the 
vacuum manifold is not simply connected (3] , or as funda- 
mental objects in string theory Q. A network of cosmic 
strings may form in thermal phase transitions 0] , at the 
end of hybrid inflation [l^, [U , or at the end of brane 
inflation when brane and anti-brane annihilate [l^ . [T^ . 
They remain of immense cosmological interest through 
their appearance in these high energy physics models, 
and further motivation is provided from the enhanced flt 
of the ACDM model to Cosmic Microwave Background 
(CMB) data when cosmic strings are included [3, 
(see also [11 [H [11 [H HI ) . Calculations of the CMB 
signal at small angular scales [U, and in the polarisa- 
tion B-mode [1^, [l^l show that future CMB observations 
will further constrain cosmological models with strings. 

Cosmic strings form networks of infinitely long string 
and loops, where a string can be called "infinite" in cos- 
mological terms if it is larger than the horizon. Infinite 
string takes the form of a random walk with correlation 
length ~ ^. The network evolves in a self-similar man- 
ner, keeping ^ at about the horizon scale, an important 
dynamical feature known as scaling. Scaling means that 
the energy density of infinite strings decreases as {t 
is cosmic time), and thus constitutes a constant fraction 
of the total. 



'Electronic address: 'm.b.hindmarsh@suss ex.ac.ukl 
tElcctronic address: s.r.stuckcy@susscx.ac .uk] 
J Electronic address: |n.bevis@imperial.ac.uk| 



Since the original string network scaling paradigm was 
introduced [3, [H, [2^ a broad picture of the cosmological 
evolution of string networks has emerged, using a mix- 
ture of calculation, numerical simulation, and analytic 
modelling. Yet, a major unsolved problem is the even- 
tual destination of the energy in the infinite strings. This 
is of notable importance and greatly limits our ability to 
constrain string scenarios via their decay products, such 
as gravitational waves or energetic particles. 

In the traditional picture, string will loop back on it- 
self and undergo a self- intersection event, which then re- 
sults in the loop being cut from the long string. Strings 
have tension, so loops tend to collapse inwards and can 
shrink to a point. In their final moments, as their radius 
becomes close to the string width, they give up their 
energy in a burst of particle emission, which in principle 
might be detectable via cosmic rays. However whilst they 
are shrinking, they would oscillate and emit gravitational 
waves. Conventionally it is taken that the loops are large 
relative to the string width and that these gravitational 
waves take away most of the energy. 

Numerical simulations using the Nambu-Goto approx- 
imation [27i, |28|, |29| , in which the strings arc modelled as 
relativistic one-dimensional entities, seemed to support 
this picture. Copious loop production was observed at 
scales a small fraction of the cosmological horizon, along 
with the presence of small-scale structure in the correla- 
tion functions of the long string. The small-scale struc- 
ture is related to the creation of small-scale loops, but 
progress on understanding the connection has been slow. 

It is crucial to establish the dominant length scale of 
loop production, as it controls both the amplitude and 
frequency of the gravitational wave signal, and the frac- 
tion of energy going into ultra-high energy cosmic rays. 
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FIG. 1: A snapshot from a simulation of a string network 
in the Abelian Higgs model. Lines show the centres of the 
strings, and the shading on the top and bottom faces repre- 
sents magnetic field energy density and Higgs field potential 
energy density respectively. 
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FIG. 2: The comoving network length scale ^ (roughly corre- 
sponding to the radius of curvature of long string) for a 1024^ 
lattice simulation in the matter era. The linear behaviour 
with conformal time r is evidence for scaling extending to 
values of ^ nearly two orders of magnitude greater than the 
string width, which is of order unity in simulation units. See 
Section [Till for more information. 



The current debate is summarised in Section V. Here 
we merely note that there is a great deal of uncertainty 
over the loop population of cosmic string networks and 
the Minkowski-space Nambu-Goto simulations have been 
used to argue that loop production might even peak on 
the scale of the string width [3^ . 

If this is true, it is necessary to go beyond the Nambu- 
Goto approximation and calculate with the underlying 
theory. The Nambu-Goto approximation also breaks 
down at kinks (discontinuities in the string tangent vec- 
tor) and cusps (points where the tangent vector van- 
ishes and the string doubles back), which are universal 
and common features of a string network. The under- 
lying theory for solitonic strings is a quantum field the- 
ory: Fortunately, quantum corrections appear to be small 
[3l| . and classical field theory should be a good approx- 
imation. Numerical simulations in the classical Abelian 
Higgs model [H, [3^ showed that infinite string does in- 
deed scale by losing energy into gauge and Higgs radia- 
tion (see Fig.[T]), although it was not established whether 
the decay proceeded via short-lived loops at the size of 
the string core width, or directly from the long strings 
themselves. In any case, no sign of copious loop produc- 
tion was found. In this paper we present evidence that 
direct radiation is much more important than small loop 
production for Abelian Higgs string networks. 

The idea of direct radiation raises a puzzle: in order to 
create radiation in a mode with mass M, the field must 
be oscillating at a frequency lo ^ M . A smooth string 
curved on the horizon scale H^^ is constructed from field 
modes with frequencies to ^ H. In view of the mismatch 
it was argued that gauge and Higgs radiation must be 
negligible [33 |. Numerical simulations of smooth strings 
|35| show that radiation is indeed suppressed by a factor 
exponential in the ratio of the curvature radius to the 
string width. Nevertheless, simulations exhibit scaling 



behaviour at times when the network len gth scale exceeds 
the string width by a factor ~ 50 [H, fssl . [36j , and for 
example in Fig. [5] we see that this is confirmed to ^ 100 
in our larger simulations. 

It should be noted that a similar puzzle was presented 
by the small size (compared with the horizon scale) 
of loops in Nambu-Goto simulations. Recent work by 
Polchinski and collaborators P, H, [sj has resolved the 
problem by showing how small-scale structure can give 
rise both to apparently smooth strings and to loop pro- 
duction at the small-scale cut-off on the string network. 

Motivated by the success of the approach by Polchinksi 
and collaborators [H, 0| , we analyse the small scale struc- 
ture of cosmic strings in the Abelian Higgs model through 
the two-point tangent vector correlator and the loop 
distribution function. We find excellent agreement for 
the two-point tangent vector correlator, showing that its 
slope at short distances and the mean square velocity are 
related as predicted by their model but with mean square 
velocities found to be significantly lower in Abelian Higgs 
simulations than in Nambu-Goto simulations. The pres- 
ence of small-scale structure on Abelian Higgs strings 
leads us to propose that the high frequencies required 
for direct radiation are generated in much the same way 
as small loops. 

However, the prediction for the form of the loop pro- 
duction function is much less successful, even after tak- 
ing into account the fact that in field theory simulations 
loops lose energy and shrink at a constant rate. From evi- 
dence found of fragmentation of horizon size loops we are 
prompted to compare our data to the model of Scherrer & 
Press 3 and this seems to provide a better explanation. 

We also investigate the speed with which the network 
relaxes to scaling, and in particular how quickly the small 
scale structure appears. We find evidence that the re- 
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laxation time goes roughly as the initial network length 
scale. Furthermore, we present evidence that small scale 
structure grows from the bottom up, rather than result- 
ing from a cascade from large scales as has hitherto been 
assumed 

Our paper is laid out as follows. A review of the tra- 
ditional picture of energy loss from cosmic strings is pro- 
vided in Sec. |TT] then the particulars of the field theory 
simulation for Abelian Higgs strings are outlined in Sec. 
mil Results of the analysis of the two-point correlation 
fimctions and loop distributions are discussed in Sec. IIVI 
and Sec. El We summarise our results in Sec. |VT] and 
present a picture of string network evolution which rec- 
onciles our results with those derived from Nambu-Goto 
simulations. 



II. ENERGY LOSS MECHANISMS 

Conventional "solitonic" cosmic strings are produced 
in a symmetry-breaking phase transition at a scale 77, in 
a field theory with coupling constant g, and typical mass 
scale m ^ grj. They have linear mass density ^ ^ rf and 
have width 5 ^ m~^. For example, strings produced at 
the Grand Unification energy scale have 77 ^ lO^^GeV, 
and gravitational coupling G/i ~ (ij/mpi)'^ ~ 10~®. 

A particularly important feature in string energy loss is 
the presence of cusps, points where the speed of the string 
becomes the speed of light. This produces a singularity 
in the shape of the string where the tangent vector to the 
loop vanishes. The string is bent back on itself and has 
the opportunity to interact with itself over the length 
of the cusp region. Also important are kinks, formed 
during intercommutation when there are sharp changes 
in the tangent vector. These discontinuities then resolve 
themselves into kinked waves travelling in both directions 
away from the intercommutation site. Numerical simula- 
tions based on the Nambu-Goto equations show a large 
number of intercommutations, forcing a build up of small 
scale structure. When oppositely-moving modes interact 
they can emit gravitational radiation, and the resulting 
back reaction [38[ can smooth the kinks, as shown in nu- 
merical simulations [s^ . 

Conventionally, three main energy loss mechanisms 
have been considered, all of which take place on loops 
which have broken off from the long string network. Per- 
turbative production of particles from the coupling to the 
Higgs field for strings with m <C /i^^^ can be calculated 
to produce a flux of high energy particles with power (40| 

Pp = i^5r\ (1) 

where I is the size of the loop. Where string is overlapping 
in cusp regions, it can annihilate releasing energy in the 
form of high energy particles. The original estimate [4l| 
was corrected in [l^l who obtained 

Pe OC /i^l/^rl/^. (2) 



The power emitted through gravitational radiation, P^, 
from a sizeable oscillating string loop of length I and mass 
m ~ ^bl can be estimated from its quadrupole moment 

P,~g(^)'^Gc.^/2ocGm^. (3) 

Much additional work has been done on the produc- 
tion of gravitational wave bursts at the sites of cusps on 
striiigs with and without additional small scale structure 
[45I . m, I43, Ell . Cusps on strings with small-scale struc- 
ture are also sources of intense loop production 0, [4^ . 

Purely based on the length dependence of these re- 
lations for the power output; gravitational radiation 
Pg is the dominant decay channel for loops of length 
I > S{Gii)~^, where S is the string core width. Below 
this length, cusp annihilation is is dominant. So it is 
of considerable importance to discover the distribution 
of different sized loops produced and evolving in the net- 
work as their size determines decay modes and the radia- 
tive by- products. This has crucial bearing on estimates 
of gravitational radiation or flux from high energy cos- 
mic rays that we may be able to detect from a network 
of cosmic strings. 

In the conventional cosmic string scenario the loop dis- 
tribution is measured in simulations using the Nambu- 
Goto approximation and neglecting gravitational back- 
reaction. This is justifiable when the radius of curva- 
ture of the string is much greater than the string core 
width so a string can be assumed infinitely thin, and 
when the weak field approximation holds, Gfi ^1. It is 
clear that in the presence of kinks and cusps, the Nambu- 
Goto approximation is strictly not justified, but it is as- 
sumed that gravitational radiation back-reaction will act 
to smooth the strings on a scale I ~ (Gy^t) ^"'"^-^t [s^l , where 
X is a small parameter defined below, and that large-scale 
properties will be reproduced correctly. 

In the absence of back-reaction, Nambu-Goto simula- 
tions show that loops are produced at a small constant 
physical scale, which is most likely the initial correlation 
length of the network [sO, [Ml, H^] . In Ref. [5l| there is a 
claim that there are signs that this scale is growing, while 
Ref. [5^ emphasises the significance of an apparently sta- 
ble population of loops with sizes I ~ O.lt, arguing that 
the peak at the initial correlation length will eventually 
disappear. 

In view of the automatic breakdown of the Nambu- 
Goto approximation when strings intercommute and pro- 
duce kinks, it would seem sensible to simulate strings 
in an underlying field theory. One resolves the kinks 
and cusps correctly, and includes classical radiation as a 
form of energy loss. While the conventional arguments 
given above emphasise gravitational radiation, omitted 
from all simulations, one should be able to see the other 
forms of energy loss and to check their scaling with loop 
size. However, as outlined in the introduction, previous 
field theory simulations [s^ . [ssj found a scaling infinite 
string network without a significant population of loops. 
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core width would quickly evaporate into classical radia- 
tion, and despite a large production rate would have a 
very low number density, which is easily estimated to be 
a few per Hubble volume [3^. In this scenario, field the- 
ory strings are behaving like Nambu-Goto strings in that 
loops are being produced at the smallest physical scale, 
which is in one case the string width, and the other the 
initial correlation length. While we see core loops, we 
would expect them to be associated with long string, and 
as we will show, we are unable to find such a correlation. 
Their number density is also too low to account for the 
energy loss from long string. It may be that energy is 
being broken off in lumps which are too small to register 
as loops at all. 



FIG. 3: Change of total comoving length of string / with 
conformal time r for late time radiation era simulations of 
box side size L — 256A2; with Aa:: = 0.5. The energy loss 
from a loop of string is linear with constant gradient until the 
loop fragments and evaporates (examples shown in black). 
Once only 2 straight strings remain in the box (3 examples 
shown in grey), the linear energy loss stops and each string 
length oscillates around the box crossing length; L, Lypl or 
L\f% depending on the direction the string wraps the box. 



There appears to be an energy-loss mechanism allowing 
the strings to scale, which has a different length depen- 
dence to any of those outlined above. 

To see this, consider a string network whose spatial 
distribution scales with the horizon size such that the 
length of string in an horizon volume is ~ The energy 
density of the string network hence varies as ps ^ m/^^, 
which then requires the network to lose energy at a rate 
Ip'sI ~ m/^^- Hence the rate of energy loss from a piece 
of string of size ~ t is Pg ~ /j,. Thus a scaling network 
requires an energy loss mechanism which is independent 
of the size of the string, like gravitational radiation, but 
a factor (G/i)~^ stronger. One implication of the mech- 
anism is that the length of a loop of string will shrink at 
a constant rate of order unity. 

We verify this by running 256"^ radiation era simula- 
tions for a very long time until there is only a single 
shrinking loop or a pair of straight strings winding around 
one or more directions in the simulation volume, which 
has periodic boundary conditions. The change in to- 
tal comoving string length with conformal time is shown 
in Fig. [3] to be linear. By inspection one can see that 
for strings that are shrinking, dl/dr ^ 0(1) as claimed. 
Once two strings in the box become straight their length 
oscillates around the box crossing distance, which is the 
behaviour expected in the Nambu-Goto approximation. 

The detailed mechanism for this strong energy loss is 
not well understood. Attention has been focused on the 
production of 'core' sized loops [H, [s^, which would 
nicely connect the Nambu-Goto and field theory simu- 
lations. Core loops produced at the size of the string 



III. STRING SIMULATION METHOD 

A. Abelian Higgs model 

As discussed in the introduction, we characterise the 
small-scale structure and loop distribution of strings via 
simulation of the Abelian Higgs model [1, d, [1]. This 
is the simplest field theory to contain local (or gauge) 
cosmic strings and has Lagrangian density: 

C = (D.q^nD^cl^) - -^F.^F^'^ - ^(|0p - r^Y , (4) 

where (/> is a complex scalar Higgs field, = lA^ is 
the gauge- covariant derivative and F^i, = dfj^A^ — d^A^ 
is the usual field strength tensor. String solutions are 
well known in this model [H^; the phase of the Higgs 
field winds around 27rn {n € Z^) as a closed loop is tra- 
versed through space and (p is forced to depart from the 
vacuum manifold over a tube of radius ~ l/\/X7]. The 
gauge field acts to compensate the winding but this re- 
sults in a pseudo-magnetic flux tube of radius ~ l/e??. 



B. String width control 

These two fixed physical length scales must be resolved 
in any simulation of the string network but, in an expand- 
ing Friedman-Robertson- Walker universe, they rapidly 
fall away from the other length scales that must also be 
resolved. Using comoving coordinates the horizon size is 
given simply by the conformal time r and is the relevant 
length scale at which the super-horizon tangle of string 
begins to straighten and decay. However, the comoving 
string width decays as the reciprocal of the cosmic scale 
factor a and therefore as in a radiation-dominated 
universe and under matter- domination. We there- 
fore require to resolve two scales which diverge as un- 
der matter domination, which would, in principle, limit 
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us to very short periods of simulation-'^. 

However in Ref. [1^1, a technique was demonstrated 
via which the couphng constants e and A were raised to 
time dependent variables: 

A = Aofl-'^i-") (5) 
e = eofl-*^-'^' (6) 

in order that the comoving string width r behaves as: 

r (X a"". (7) 

That is, K = 1 gives the true dynamics while k = yields 
a comoving string width, which is particularly convenient 
for simulation. The dynamical equations derived upon 
variation of the action were then: 

4> + 2U- D,D,<t> = -a^'-^m' - V^)^ (8) 
Foj + 2(1 - k)^Foj - = -2a2«e2XTO[0*i^,0] (9) 

(in the gauge Aq = 0) with the Aq variation yielding the 
quasi-Gauss' Law constraint: 

- d,Fo., = ~2a^^epm[cj)*<j)]. (10) 

Although these field equations do not conserve energy if 
K 7^ 1 (because the action is not time- translation invari- 
ant), the dynamics were shown in Ref. (36j to be insen- 
sitive to K. Indeed the difference in their results for the 
two-point correlation functions of the energy-momentum 
tensor between k = 1 and k = in the radiation era, 
where k = 1 was achievable, were found to be slight 
while their results for the CMB power spectra, which are 
dominated by the matter era, were found to be similarly 
insensitive to k over the range to 0.3, with the later be- 
ing the largest value at which reliable CMB results could 
be obtained. 

Here we use the equations of motion Eqs. ([8]) and ([9|) 
with K = throughout, although we also check our re- 
sults using simulations with k = 1 for the radiation era 
(only), and find changes that are negligible. Note, how- 
ever, that our equations of motion at k = are not pre- 
cisely the same as those used in Ref. [s^ . 



LAH, which is an extension from Minkowski space-time 
to flat FRW universes of the standard approach of Ref. 
[sl ]. This preserves the Gauss' law constraint to ma- 
chine precision. Initial conditions were chosen following 
Ref. [36[ in order for a string network to emerge rapidly 
and the Gauss constraint obeyed. To achieve the latter, 
we set to zero the gauge field and the time derivatives of 
both the gauge and Higgs fields. To achieve the former, 
we set the simulation start time such that the horizon 
is comparable to the (uniform) lattice spacing Ax and 
therefore the phase of scalar field is an independent ran- 
dom variable on each lattice site, while we set \(j)\ = rj. 
We employ periodic boundary conditions and therefore 
the fields can be evolved forward reliably until the half- 
box crossing time for light. 

We use a lattice spacing of Ax = 0.5 and set 77 = 1, 
Aq = 2 and eg = 1 , which together guarantee that strings 
are resolved for all a when k = and for a < 1 in our 
radiation era k = 1 check. Note that the above scalings 
leave the ratio A/2e^ constant and we therefore study 
the model at the Bogomol'nyi value [5^, which yields 
equal scalar and gauge masses. The simulations were 
performed using the UK National Cosmology supercom- 
puter [5^ , with parallel processing enabled via the LAT- 
field library [s^l , using a lattice size of 512^ and averaging 
over 20 realisations. Additional simulations using a 768"^ 
lattice were performed on the University of Sussex HPC 
Archimedes cluster. 

We have been able to access a dynamic range of simi- 
lar order to Nambu-Goto simulations performed by other 
groups working on small scale structure issues but due to 
the differences between the simulations the measures for 
dynamic range that are suitable in one case are ambigu- 
ous in another, making this a difficult comparison. On 
lattice volumes with Ax = 0.5, a network of strings of 
width (5 ^ 1 is fully formed from conformal time ~ 20 
and is simulated until ry = NAx/2, when the periodic 
boundary conditions of the simulation volume can poten- 
tially be felt. Checks of the full energy- momentum tensor 
indicate that scaling is achieved at around Tsc ^ 64. 



C. Simulation Specifics 

Eqs. ^ and ([9|) are discrctised onto a lattice using 
the procedure described in Ref. js^] and hereafter called 



^ That our simulations resolve the string width limits them in size 
to being far smaller than the horizon size at matter-radiation 
equality and therefore our simulations are also limited to much 
earlier times. However, we can simulate a matter dominated 
universe at very early times and use scaling to make statements 
about a matter dominated era at later times. 



One measure of dynamic range is ^(t/)/(5 ~ 50, which 
contains measurable quantities in our simulation. This 
can be compared with the ratio of the initial and fi- 
nal correlation lengths, (~ 100 [13, [131) in Nambu-Goto 
simulations although the initial correlation length has 
no straightforward meaning in our simulation. Another 
measure is the ratio of the final time to the time at which 
the network achieves scaling. For us, t^c ~ 64 gives us 
a dynamic range of about 2 (in conformal time) for our 
512^^ simulations and about 3 for 768'^ lattices. Nambu- 
Goto simulations [13| estimate a dynamical range of 5 
from the scaling of the energy density of long string in 
the radiation era. 
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D. String length measurements 

At intervals during the evolution we record the coor- 
dinates of lattice plaquettes around which the phase of 
(j> has a net winding number^. As a first approximation 
we then take it that a segment of string having length 
Ax threads through each plaquette of 2tt winding and 
joins the centres of the lattice cells on either side. From 
these segments we can then construct the path of the 
string, although since it is composed of an array of per- 
pendicular lengths the string length is overestimated. We 
do not attempt to smooth the paths in order to com- 
pensate, as in Refs. (s^ . Issl . ISa ]. but instead apply a 
Scherrer-Vilenkin correction of tt/Q [g^I to such length 
measurements. This factor is derived from the length of 
the line (/> = for a two-component Gaussian random 
field (j) and so will not be completely accurate for our 
dynamic string network. However, the results for the av- 
erage string length density are in approximate agreement 
with that measured using the Lagrangian density Fig. U) 
This second method makes use of the fact that pertur- 
bative radiation makes zero contribution to C, while a 
static straight string contributes at — /i per unit length 
density. Since £ is a four-scalar and length density picks 
up a 7-factor upon a Lorentz transform, then — J C(Px/ fi 
is a measure of the (invariant) string length. 

Rather than use the (comoving) string length density 
directly, we instead compare with the network length 
scale ^, defined as: 




where V is the reference volume and L the string length 
within it. For a scaling network ^ oc r. 

In Fig. m we plot both the Lagrangian measure re- 
sult and the winding measure result ^ (with no sub- 
script since it is our default measure), which reveal a 
linear behaviour ^ oc (r — t^=o) after an initial transient, 
consistent with the expected scaling. As pointed out in 
Rcf. [3^, there is nothing fundamental about the value 
of T^=o a-nd it is simply an artefact of the initial condi- 
tions. Indeed, this can be seen in Fig. [5l which shows 
additional results from two runs in which an artificial 
damping phase (similar to that used in Ref. [6lj) was 
employed for an extended time. As can be seen in the 
figure, when the damping is released, both these runs 
show ^ rapidly reverting to approximately the same gra- 
dient as the undamped run. Reference [36| also observed 
scaling with r — tj=o in the two-point correlation func- 
tions of the energy-momentum tensor, so there is no evi- 
dence from the simulations that this scaling is a transient. 



^ While a winding of the phase is gauge-invariant in the contin- 
uum, on a lattice it can be removed by a finite gauge transform. 
Therefore we employ the gauge- invariant measure of Ref. ISSlI 
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FIG. 4: Network length scale ^ in a radiation era simulation 
calculated from Eq. is in solid dark grey with fit over 
scaling times (r £ [64, 128]) giving ^ oc 0.31r (the dotted line 
shows ^ before rescaling total string length using the Scherrer- 
Vilenkin factor). 5£ calculated directly from the Lagrangian 
is in solid light grey with fit oc 0.29r. 
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FIG. 5: Network length scale ^, showing linear behaviour 
with conformal time once the network settles into a scaling 
regime and memory of initial conditions is lost. ^ (black) and 
(grsy) are shown for radiation era simulations according 
to Eqs. (|8} and (|9| with parameters as stated in Sec. IIII CI 
Curves for smaller ^ from simulations with a constant damp- 
ing term which is not switched off until t — 75 and r — 100 
respectively showing the speed with which the network length 
scale resumes the same scaling evolution with ^ oc 0.31r. 



These energy-momentum correlators show that the net- 
work demonstrates scaling in a 512^ simulation over the 
conformal time range 64 < r < 128, which will be re- 
ferred to as the scaling epoch. 
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IV. CORRELATION FUNCTIONS 

A. Tangent Vector Correlators 

A 2-point correlation function from LAH simulation 
data is compared to analytic results which are reviewed 
in a manifest form with a view to introducing the vari- 
ables and functions under discussion. For a thorough 
analysis refer to [l|, 0, [s^l ■ The equations of motion are 
reformulated in terms of the left and right moving unit 
vectors p and q defined in terms of the position vector 
x((T, t) where a is the string coordinate. 



p((T, r) = X + ix' and q(cr, r) = x — ^x' 



where dots are derivatives with respect to r, primes with 
respect to a and 



(13) 



1-^2 



thus giving the equations of motion in the transverse 
gauge m 



P 

q 



= -f(q-(p-q)p) 
= -i(p-(p-q)q)- 



(14) 
(15) 



From the coordinate locations of the string extracted 
from LAH it is simple to compare the Euclidean distance 
between 2 points on the string and the separation along 
the string coordinate. 

The longest string at a set of equally spaced times 
throughout the scaling epoch in the simulation is isolated 
for analysis in both the radiation and matter dominated 
eras. The comoving distance along the string s = J eda 
along a string coordinate length ct = cti — (T2 is compared 
to the Euclidean distance, r, between ai and tT2- Along 
each loop of string, the mean square Euclidean distance is 
averaged over many starting points a few string segments 
apart around the loop, thus creating a 2-point function, 

{r'i<y,r))^ r rdV(x'(ai)-x'(a2)). (16) 
Jo Jo 

For consistency with notation used elsewhere P, [5l| 
and using x' = § (p — q) in terms of the definitions of Eq. 
(fT^ we denote 



corrx (cr,T) 



(x'(ai)-x-(a2)) 
(x'(0) • x'(0)) 

(P(o-l) ■ P(g2)) - (P(o-l) ■ q(0-2)) 

2(1 -(x^)) 



(17) 



Correlations for opposite moving modes (p(cri) • q(cr2)) 
will be subdominant to the build up of those moving 
in the same direction along the string (p(cri) • p(ct2))- 
Identifying that — p(cr) • qlc) = 1 — 2x^, averaging over 
an ensemble of segments and many bubble times allows 



an approximation to first order of — (p(cri) • q(cr2)) = 
1 — 2(x^) for small s giving 

2(1 - (i2))(l - corrx(a,T)) = 1 - (p(a) • p(0)). 

The equations of motion can be rewritten in terms of the 
evolution of (p((t) • p(0)). 



2d , 



a.(l-p(a).p(0)) = (l-2(x^))(l-p(a).p(0)) (18) 

a 

which integrates to the form 

(l-p(a).p(0))=/(a)r-2-(i-2(x^» (19) 

with v defined as scale factor a (x t'^ . The correlator must 
be a function of s /t if it is to exhibit scaling. Given that 
s ~ ecr, we need the time dependence of of e (Eq. [13)) 



- = -2V^eocr-2'^<^'>. 



Hence 



Given that Eq. has a power law form 
(1-P(a)-P(0))(x (^)''^ 



we find 



2X 



(20) 



The tangent vector correlator should also scale 



1 - corrx(cr,T) = 



s\2x 



(21) 

x'.x')(i-A(0'') (22) 



In order to test the prediction, the second derivative of 
the 2-point function (r^(s,r)) is taken numerically to 
find a fit for the parameters in this analytic result. A 
least squares fit is used to optimise the parameters in 
the function for 1 — corrx taking a nominal standard de- 
viation equivalent to the length s to weight smaller s 
appropriately on the logarithmic scale. Noise is reduced 
by averaging the result over 20 simulations and addition- 
ally averaging over a block of As = 20 and centring the 
mean value. It should be noted that the smoothing pro- 
cess was not found to alter the parameters or the shape of 
the 1 — corrx function, notably at small s, but allowed the 
least squares fit to converge more quickly. These func- 
tions are shown for the radiation era Fig.lHland matter era 
Fig. [7] simulations. The correlations over long distances 
refiect smoothness on large scales and a clear power law 
demonstrated on smaller scales from a few times the cor- 
relation length down to string width scales. The average 
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V 

0.5 

1 

1.5 
2 



(x-x'} A 2x 

1.0 0.86 0.44 

0.97 0.40 0.62 

0.96 0.85 0.71 

0.92 0.30 0.90 



{X^)2x 

0.23 
0.26 
0.27 
0.27 



(x!)g 
0.22 
0.27 
0.27 
0.26 



TABLE I: Table to show parameters found from 2 point cor- 
relation data fitted to predicted power law Eq. ([21]) of Ref.[H. 
Correlation functions plotted in Figs. [5] and [7] The implied 
mean square network velocity, {x^)2x is calculated from the 
exponent according to Eq. (|20|l . These values are compared 
with the direct calculation of network velocities in Fig. |S] and 
quoted for (x^)g, calculated from the field gradients. 




-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 
log(s/5) 



FIG. 6: Solid black lines indicate 2 point tangent vector corre- 
lation functions in the radiation era for 7 equally spaced times 
throughout the scaling epoch (r £ [64, 128]) with dashed grey 
fit calculated from an average of all 7 sets of parameters fitted 
up to log(s/5) = 0.5. 
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FIG. 7; Correlations and fit for the matter era as in Fig. [6] 
with fit performed on curve to log(s/^) = 0.3. 
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FIG. 8: From bottom to top (x^) is calculated for v = 
0.5,1,1.5,2. Crosses mark values calculated from 2x expo- 
nent according to Eq. (|20|) : dark grey shows (x^)f calculated 
from Eq. (|23p : light grey shows {x^)g calculated from Eq. 



parameter values over 7 seven equally spaced times in 
the scaling epoch of the simulation have also been cal- 
culated for v = 0.5 and v = 1.5 as an investigation of 
transitions between eras. Results are given in Table HI 
The parameter 2x provides an estimate for the average 
network velocity squared (x^) via Eq. pO|) and its value 
can be cross-checked with that calculated directly from 
the simulation Fig. [51 Two methods are used for estimat- 
ing the average network velocity. One uses the electric E 
and magnetic B components from the field strength ten- 
sor and the other the canonical momentum and spatial 
gradients of the field, 11 = </) and D0 = d(j>/dx + iA(f>. 
Denoting the estimators (x^)f and (x^)g, we have 



and 



E| 



(D<; 



(23) 



(24) 



where jp and are calculated using (x^)f and 
respectively, and the subscript C denotes a Lagrangian 
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weighting of a field X according to 

We see that the measured velocities agree well with those 
inferred from the slope of the correlation function, strong 
evidence that the model of [H describes the dynamics of 
long string in the Abelian Higgs model well. 

Finally we note that our root mean square (RMS) ve- 
locities are approximately 0.5 in both the matter and ra- 
diation era. These are significantly lower than measured 
in Nambu-Goto simulations, about 0.66 in the radiation 
era and 0.61 in the matter era. This is likely to be a re- 
sult of backreaction from massive radiation, not included 
in Nambu-Goto simulations. Our velocities are, however, 
in agreement with the uncorrected RMS velocities of Ref. 
[3^, which are measured in a field theory simulation us- 
ing the positions of the zeros of the field. 



B. Initial Conditions and Relaxation to Scaling 

It is important to test for any dependence of our re- 
sults upon the initial conditions chosen and to fully ex- 
plore the approach to scaling. To achieve these two goals 
we have performed additional simulations with an initial 
phase in which the Hubble damping term 2d/a in Eqs. 
dH) and ^ is replaced by a constant, which we set to 
unity. This initial phase of artificial damping continues 
until a time when the Hubble damping would have been 
far smaller and hence string velocities are heavily reduced 
by the time we switch over to normal Hubble damping 
to complete the simulation. This gives us an alternative 
initial condition in which the string network is smooth 
and slowly moving while radiation is negligible. 

The effects on the network length scale ^ for simula- 
tions on 512'^ lattices, seen in Fig. [SJ shows the rate of 
growth of ^ is heavily retarded during the initial phase 
with ^ cx r^/'^ as expected for over-damped motion [63 |. 
Once the constant damping is switched to Hubble damp- 
ing, the evolution of ^ quickly transitions to the same 
growth rate seen in the primary simulations. 

The effects on the correlation functions are shown in 
Figs. [QlandflOl For clarity, the log of 1 — corrx is plotted 
against the log of the comoving length along the string s, 
which separates the curves so the changes in gradient are 
more apparent. It is clear from the tangent vector cor- 
relations that when the evolution is artificially damped, 
the small scale structure is quite different. The initial 
slope for log(s) < 1.4 is close to unity, as expected for 
smooth strings, which is maintained for a short time. 
Then 1 — corrx begins a transition towards a new expo- 
nent as scaling is attained for the new undamped (Hubble 
damping only) evolution. The power law changes first at 
small scales and moves up to larger scales as the evolu- 
tion adjusts to the new conditions. Once the exponents 
have relaxed, the small scale structure behaves and scales 
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0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 2.2 2.4 

log„(s) 

FIG. 9: Log of the tangent vector correlation function 1 — 
corrx against log of the comoving string coordinate length s 
for a radiation era simulation with constant damping until 
r = 55 on a 768^ lattice. Black curves show the correlation 
function at times before, at and immediately after damping 
is released to ordinary Hubble damping. In dark grey are 
correlators at times when then the exponent is in transition 
from the constantly damped evolution but before scaling is 
achieved (r = 62, 74, 80, 88). Scaling is re-established for the 
correlators plotted in light grey for times r > 96. 
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log(s) 

FIG. 10: Radiation era simulation on 768^ lattice with con- 
stant damping until r = 100. The black curve shows 
log(l — corrx) at r = 96; before artificial damping is 
switched off. In dark grey are correlators shown at times 
r = 104,112,120,128,136 when then the exponent is in 
transition. Scaling correlators are shown in light grey for 
r = 144, 160, 176, 192. Noise is greater for these curves com- 
pared with Figs. |6] and [7] as data is taken from just one sim- 
ulation. 



just as in an undamped simulation with no recollection 
of the initial damping conditions. 

The exponents are found by setting (x' • x') = 0.97 as 
found for the radiation era simulations on a 512'^ lattice 
which are also consistent with that found for a single 768'^ 
simulation. So in this case only a 2 parameter fit is made 



10 



- noDamp 

- Damp40 
Damp55 
DamplOO 
Damp200 
Damp300 




FIG. 11: Relaxation of the parameter 2x after constant damp- 
ing is switched off in a number of 768'^ radiation era simula- 
tions with Toff = 0, 40, 55, 100, 200, 300. 



FIG. 12: Grey hues show correlations for k — and black 
for K = 1. In the top figure the comparison is made at time 
T — 85 and at the bottom is the comparison at time r = 115. 



for Eq. ([22|l : 2x and A. A is found to be constant and 
consistent in both a damped and undamped evolution 
so we discuss the relaxation time scale in terms of the 
time taken for the gradient parameter, 2%, to transition 
back into scaling. There is little effect on the value of 2x 
once the new scaling value is reached after damping and 
becomes approximately 0.6, the same as found during 
the scaling epoch for a completely undamped simulation. 
The profile of this relaxation is compared for different 
damped simulations and shown in Fig. 1111 

Although not conclusive, there is evidence that the re- 
laxation time, from the moment damping is turned off 
Toff, is approximately ^. One can see the trend most 
clearly in the value of t — Tqs at 2x — 0.7. This is expli- 
cable if the transition is triggered by intercommutation 
and the generation of kinks, which would have the ob- 
served effect of reducing the correlation at small scales. 

Figs. [51 and [TUl show the correlators settling into a scal- 
ing regime behaving much as the evolution with Hubble 
damping only. We recall also that the network length 
scale ^ returns to the expected linear behaviour. Fig. [51 
over perhaps an even shorter time frame. We interpret 
these features as good evidence that calculations can be 
trusted to not depend on initial conditions. 



C. Comoving String Core Width 

Simulations conducted in the radiation era are possible 
with the true equations of motion using k = 1 in Eqs. 
([5]) and correlation results for an average of 5 simulations 
for the correlation function results are tested against n = 
results. Fig. [T^ shows a comparison of the correlations 
at two different times in the simulation for both the n = 
and the k = 1 cases. The difference in the results 
is surprisingly insignificant for this calculation and no 
correction is felt necessary. 



D. Testing Gaussianity 

The 4-point correlation functions (r^(s,T))^ and 
((r^(s,T))^) are also calculated to test for gaussianity. 

Denoting each spatial component of r as = (cti ) — 
Xi(cr2), with CTi some initial base point for the measure- 
ment, the fourth moments formula is written 

= d^jSjk{{A^Aj){AkAi) 
+{A,Ak){A^Ai) 
+(A,Ai)(AjAfc)) (25) 

Then gaussianity would allow contraction on all pairs 
from Wick's Theorem so that ({A.,Aj)) = ^5ij{A'^) and 
the ratio of the 4 point functions should behave as 

Fig. [T3] shows for the radiation era that on all scales 
the ratio of the 4 point correlations is not constant and 
not 5/3. The 4-point correlators in the matter era (not 
plotted) behave in a similar way. 

In the model for loop production proposed by Polchin- 
ski and collaborators P, [3, [131, argued that non- 
gaussianity does not affect the power laws derived for 
the 2-point correlation function or the loop production 
function. We have confirmed that the 2-point correla- 
tion function is in accord with their model, but in the 
next section we will see that the loop production func- 
tion is not. If we accept the arguments of Polchinski et 
al, which are based on scaling, another reason must be 
found to account for the difference. 
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FIG. 13: The ratio of the 4-point correlation functions to the 
Gaussian prediction in the radiation dominated era is shown 
against the dimensionless distance parameter s/^. A Gaus- 
sian 4-point correlator should give a value 1: Gaussianity is 
therefore not evident at any scale. 



V. LOOP DISTRIBUTIONS 

Loop production in cosmic string networks is still a 
subject of some debate despite a number of recent numer- 
ical investigations using the Nambu-Goto approximation 
(sol . [Fll . [5^ , taking advantage of improvements in compu- 
tational facilities and algorithms, and focusing on small 
scale structure and loop production rates. 

The crucial quantities in question are the loop (length) 
distribution function and the loop production function. 
Unfortunately, the different groups measure different 
quantities, and emphasise different features, so the re- 
sults are difficult to compare. 

Those that measure the loop production function [5l|, 
[52! find that it peaks at a small scale, with a power- 
law rise (s^ . and a less prominent feature at about a 
tenth of the horizon length t. The identity of the small 
scale is not clear, but on inspection of the data [HI, [111, it 
appears to be related to (and at least no greater than) the 
initial comoving correlation length. Full scaling requires 
that the only scale in the distribution and production 
functions should be t: the peak therefore does not scale. 
Furthermore, it is found that the amplitude of the power 
law does not scale either [52j . 

Measurements of the loop distribution function on the 
other hand [s^l, show a peak at the initial numerical 
cut-off, and scaling at intermediate scales. The peak is 
understandable as a transient from the initial evolution, 
but as the distribution function is essentially the time in- 
tegral of the production function, the intermediate range 
scaling is a puzzle. 

It has been suggested that the non-scaling of the loop 
production function is a transient effect [sTI . Is^ , and that 
the peak should eventually disappear altogether [H^l or 
start scaling if only a large enough simulation could be 



performed [5l| . However, the evidence that a power- 
law with a small-scale cutoff is a real feature of Nambu- 
Goto string networks has been strengthened thanks to 
the agreement with Polchinski and collaborators' model 
of loop production P, 0, [s^l • There is also no evidence 
for scaling of the peak in the loop production function 
from visual inspection of the graphs in Refs. [Fll. Is^. 

Accepting the evidence of a power-law form for the 
loop production function, a small-scale cut-off is required 
to keep the total energy loss finite. The conventional 
string scenario demands full scaling, and invokes gravi- 
tational radiation reaction to change the loop production 
scale to a constant fraction of the horizon size, which is 
{G^Y'^^^t, according to Ref. [13 • However, there are 
no network simulations including gravitational radiation 
reaction so this is still a conjecture. It could equally 
well be that loop production really does not scale as the 
Nambu-Goto simulations suggest; this does not prevent 
the energy density of the long string network from scal- 
ing. Furthermore, if the small-scale cut-off is the string 
width [sot 7 it is necessary to perform field theory simula- 
tions in order to include the true small-scale physics. 

Previous field theory simulations [H, have not 
studied the loop distributions in any detail, but it is al- 
ready clear that their properties are very different from 
the Nambu-Goto versions. The number of loops in the 
simulation volume is substantially less, which prompted 
the suggestion [s^l that the network could lose energy to 
classical radiation directly rather than via the produc- 
tion and eventual decay of loops. Arguing in favour of 
loop production, it was pointed out in [33| that even if all 
the energy is lost to "core" or "proto" -loops (loops whose 
length is of order the string width) that the number den- 
sity would be very low anyway, approximately t~^. It was 
also conjectured that these core loops would eventually 
grow if a large enough simulation could be performed. 

In the first part of this section we test the hypothesis 
that a substantial fraction of the energy loss from long 
strings is in the form of core loops. The impression given 
by visualisations such as Fig. [1] is that direct radiation 
appears to be very important, although it is very difficult 
to tell the difference between a large amplitude excursion 
by the Higgs field and a core loop. However, if energy 
loss into core loops is important we would expect to find 
core loops near long strings, and our first test is to look 
for these correlations. 

We define a core loop as a loop with the minimum 
number of lattice segments to create a closed loop (i.e. 
4 linked segments, located as described in Sec. IIII Dl) . 
With Ax ~ 0.5, the core loop length is approximately 
the string width. We then measure the distance from 
core loops to the closest point on a neighbouring piece of 
string and the length of the string to which it is closest. 
The results are shown in Fig. [T31 Interestingly, it is seen 
that core loops lie close to other very small loops and 
that these clusters or isolated core loops lie at distances 
of order half the correlation length from long string; be- 
tween two long strings. As loops collapse they appear to 
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FIG. 14: Distance of core loops to other string in the radi- 
ation era in the top plot and matter era in the lower plot. 
Small loops lie close together in clusters in the voids between 
long strings at about half the average interstring distance, 
~ ^/2. Emission of core loops directly from long string is not 
detected. 
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FIG. 15: Number of loops per horizon volume in the radiation 
era over the scaling epoch r G [64, 128] with data smoothed 
by averaging results over blocks of lOr. Dark grey represents 
proportion by number of 'proto' loops of size the order of the 
string core width, 5; mid-grey represents loops up to the size 
of the correlation length, ^ at that time in the simulation; 
light grey represents all string that is greater in length than 
5 and considered to be infinite string in this case. 



fragment into clusters of very small loops but the lack of 
core loops close to long strings argues against small loop 
production being a significant source of energy loss from 
long strings. 

We can also check the hypothesis by looking at the 
number distribution of loops in field theory simulations, 
using a large number of runs. The length scales of interest 
are the string width (core loops) and the network corre- 
lation length, ^, defined as Eq. (fTT|) . Fig.fTSlshows cumu- 
latively the number densities of loops per horizon in the 
radiation era over the conformal time range 64 < r < 128 
when the network is scaling. The loops are divided into 
those of length 4 links (core loops), those up to length 
^, and those beyond. Core loops are seen to occupy a 
constant small fraction of the number density. In each 
of these classes the number density of loops per horizon 
appears to be constant. The number of core loops per 
horizon is very small, of order 0.1. 

We can estimate whether this is consistent with core 
loops being a significant channel of energy loss for long 
strings. If a network with comoving length scale ^ = 
decays into loops of size I, then their lifetime should also 
be /, given the shrinking mechanism outlined in Section 
im The number of loops per horizon volume should 
then be ~ l/x^. Given that ^ 0.1, there are roughly 
100 times too few core loops if they were to take a sig- 
nificant amount of energy away from long strings. 

Our result seems to be in contradiction to Ref. [s^, 
who use a fit to the Velocity-dependent One-Scale model 
[g^ to argue that loop production is significant in their 
field theory simulations. However, they do not give ab- 
solute values of the loop distribution function and so it 
is not possible to compare the results directly. One pos- 



sible resolution, explored in more detail below, is that 
horizon-size loops with lifetime ^ t arc carrying away an 
appreciable fraction of the energy. 

To study the loop distribution and production func- 
tions in more detail, we must model both the production 
and shrinking of loops. We denote the loop distribution 
function in terms of the cosmic time t and physical length 
Ip as np{t, Ip), where Ip = a(t)l and I is the comoving loop 
length, which is given in terms of the string variables e 
and (7 by Z = J edcr. Wc denote the comoving loop dis- 
tribution function in conformal time as n(r, /). Then the 
number density of loops np(t,lp)dlp in physical length 
interval [lp,lp + dip] is related to the comoving number 
density of loops n(T, l)Al in interval [I, I + dl] by 

n{T^l)Al ~ a^np{t,lp)dlp 
n{T,l) =a^np{t,lp) (26) 

The equation governing the loop distribution function 

is 

where we introduce the loop production function in phys- 
ical units Pp. We also take into account energy loss from 
loops which we shall assume takes place at a constant 
rate such that ^ = —A. We estimate A ^ 0{1) from 
the properties of the energy loss mechanism outlined in 
Sec. inP . Using Eq. ([26|) we can relate the comoving num- 



^ Gravitational radiation would give A ~ FG/i were it to be included 
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ber density distribution and the loop production function 
in comoving units 



dn 



a , dn 
a oi 



P{t,1) 



(28) 



where P(t, I) ~ a^Pp{t, Ip). 

Assuming scahng, the comoving loop production func- 
tion and number density distribution behave as Q 

n(r,/) - ^N{x) and P{t,1) - ^/(x) 

for functions N and / of the dimensionless ratio of loop 
length to horizon size x = 1/t. Rewriting Eq. ([^5]) in 
terms of N and / one obtains (with a ozt'^), 

{x + \)N'{x) + {v + A)N{x) = -fix) 

with solution 

/"OO 

N{x) = {x + A)-(''+4) / f{x'){x' + XY+^ dx'. (29) 



Numerical simulations suggest a power law for loop pro- 
duction, / cx below X ^ 1. If radiative effects can be 
neglected (x 3> A), and making the reasonable assump- 
tion that / vanishes for x S> 1, we have from Eq. 



cx/cxx". 



(30) 



For length scales where radiative effects are strong (x ^ 
A) 



N i 



,,01+1 



(31) 



To make our measurement we define the comoving loop 
number density in a length interval Al = I 



An 



n{T, l')dl' 



Figs. [TBI and [T7l show an estimate for 

taken from the average of 20 runs in radiation and mat- 
ter eras respectively. The solid black line in Fig. [16] the 
the initial loop distribution function which shows good 
a^eement with the expected power law of slope —5/2, 

For a network that has reached scaling, the analytic 
model of Ref . [l| , further refined by Dubath et al in [2| , 
proposes that loop density is dominated by recently pro- 
duced loops. They derive a function for loop production 
which they extrapolate to loop number density distribu- 
tion under the assumption that radiative effects can be 
ignored and obtain 



I \ 2x-3 



(32) 




FIG. 16: The loop number density A'^ = r'^An/Al, per co- 
moving logarithmic bin length Al is shown for 7 equally 
spaced times throughout the scaling epoch of the radiation 
era against the dimensionless ratio x = l/^. This is compared 
with the very early time r = 10 case in solid black which 
is compatible for I ^ t with slope -5/2 (shown dotted) as 
predicted in [63]. The slope predicted by the model of Ref. 
01 including radiative effects at small scales as in Eq. (|31|l is 
shown with a+1 = —1.4 (dashed black), where a is the slope 
of the loop production function. 




FIG. 17: Figure showing number density of loops during the 
scaling epoch in the matter era with slope a + 1 = —1.1 
predicted by the model of Ref. including radiative effects 
at small scales as in Eq. (|31|l shown in dashed black. 



with X defined in Eq. (PO)) . If radiative effects are con- 
sidered the exponent for the number density distribution 
function will be higher by -1-1 at small length scales by the 
arguments from Eqs. (|29p and pT|) . This is also consis- 
tent with the analytic findings of Rocha in Ref. [6^1 where 
the loop distribution model is enhanced by including of 
the effects of gravitational radiation. From the values of 
2x calculated from our simulation results and given in 
Table m this would give exponents of aR -1- 1 ^ —1.4 in 



the radiation era and aM + 1 —1.1 in the matter era 
which are compared to our data in Figs. [TOl and [T71 The 
agreement remains less than convincing, particularly for 
the radiation era. 

Olum et al [s^l calculate the loop production function 
/(.t) from their Nambu-Goto simulations. The function 
drops from a small-scale peak with a power law consistent 
with the model proposed by Dubath ct al 0: / oc x^^~^ . 
The exponent is calculated using Eg. l[20|) with veloci- 
ties taken from simulations of Ref. [3^; vji = 0.63 and 
wjv/ = 0.57. No fit values for the gradient of log/ are 
quoted in Ref. [s^ but pictures showing average gradi- 
ents of their loop production function in the matter and 
radiation era are used in Ref. to demonstrate their 
model. Exponents are hsted in Table HI 

Number densities can also be compared with Nambu- 
Goto simulations of Ref. [13] who quote a length distri- 
bution, 

xN{x) oc (33) 

They find a consistent power law over the whole range of 
their Nambu-Goto simulation with exponents p = —1.6 
for the radiation dominated era and p = —1.4 in the 
matter era. Given that there is no decay mechanism 
in Nambu-Goto simulations we can infer slopes for the 
loop production function of a = —2.6 and a = —2.4 for 
radiation and matter eras respectively, in good agreement 
with the values predicted by the model [3] for Nambu- 
Goto strings. However, even if radiative effects are taken 
into consideration these results are steeper than those 
predicted for our Abelian Higgs strings, (see Table [III . 
which are already too steep to fit our data. 

The flatter power law for loop distributions found in 
the field theory examination Figs. [TBI and fTTl can be better 
explained by invoking loop production and fragmenta- 
tion at horizon scales combined with energy loss at small 
scales. From Ref. 0] it is predicted that unusual power 
laws for the production of loops can be explained by loop 
fragmentation probabilities, q, with / oc l^^"^. For small 
loops where radiative effects are strong, number density 
distribution functions will take the form N oc x~'^'^~^^. 
Then the steepest slope possible would be —1 when the 
fragmentation probability is of order unity. 

We also see evidence of loop production at the horizon 
scale in the small feature visible in the loop distribution 
function at logx ~ 0.5, which has some similarity to the 
Nambu-Goto simulations of [5^. It is straightforward to 
check that the production of one such loop per horizon 
volume per Hubble time is sufficient to remove a signif- 
icant fraction of the energy in long strings. Given that 
these loops are losing length at a rate of order 1 as well 
as fragmenting, this is consistent with our observation of 
order one loop per horizon volume at any time. The cor- 
relation of core loops with other small loops shown Fig. 
[14] can also be explained by loop fragmentation. 

Finally, we note that the large x = 1/t behaviour in 
our number density analysis remains a puzzle, as it de- 
parts from the -2.5 slope expected outside the horizon. 
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an aM 


Abelian Higgs Prediction 
Nambu-Goto Prediction 
Nambu-Goto Measurement 


-2.4 -2.1 
-2.8 -2.5 
-2.6 -2.4 



TABLE II: Table to show comparison of simulation results for 
the exponent a obtained for a power law model for loop pro- 
duction as Eq. (|32|l in both radiation (R) and matter (M) eras. 
The first line shows our values calculated using the velocities 
{x^)2x from Table [Ij The second line shows the exponent 
predictions calculated using velocities obtained from Ref. [5l[l 
which fits well to the Olum et al Nambu-Goto simulation [53 |. 
The last line shows exponents derived from measured length 
distributions in Nambu-Goto simulations by Ringeval et al 

M- 

It maybe that the loop distribution at these much longer 
length scales is quite sensitive to the finite volume of the 
simulation [66l. 



VI. CONCLUSIONS AND REMARKS 

In this paper we have examined small scale struc- 
ture, energy loss and loop production in cosmic strings 
found in the Abelian Higgs model (AH), using numeri- 
cal simulations, and compared the results to those from 
Nambu-Goto (NG) simulations in the literature. We have 
also investigated the energy loss mechanism which causes 
Abelian Higgs strings to scale without gravitational ra- 
diation. 

We first studied the two point tangent vector correla- 
tion function (Eq. I17p , for which there is a clear predic- 
tion based on the NG equations by Polchinski and Rocha 
(PR) [H. Their model predicts that the correlator scales, 
i.e. that it is a function of the ratio of distance to the hori- 
zon scale, and that this function has a power law form 
with an exponent determined by the mean square veloc- 
ity of the string. The correlator is seen to scale very well, 
and good agreement with the power law is found right 
down to the string core width scale. The mean square 
velocity is about 0.26 in both the radiation and matter 
eras, markedly lower than in NG simulations (0.44 and 
0.37 and respectively [H, [5l|). This is likely to be due 
to back-reaction from the massive radiation. 

We test the sensitivity of the tangent vector correla- 
tor to the initial conditions by damping the strings for 
an initial interval of various duration. We see that the 
network regains its scaling form after a time of approxi- 
mately ^. In this relaxation period we see evidence that 
the build up of small scale structure appears to move 
from small scales to large. It is plausible that the trig- 
ger for the production of small scale structure and the 
relaxation of scaling is intercommutation of long strings. 
The fact that the structure appears bottom-up argues 
against small-scale structure being due to residual power 
from the Brownian random walk at super-horizon scales, 
as argued in [l|. 
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The small scale structure seen in the tangent vector 
correlations can help to explain the scale separation prob- 
lem described in Sec. which is to understand how mas- 
sive radiation of frequency around about the string width 
is produced from fields which are apparently changing 
with a frequency of about the Hubble rate. The solu- 
tion is that if strings arc not smooth on small scales, 
the radius of curvature at these scales is much smaller 
and frequency of modes of oscillations of the string much 
higher than the Hubble rate. 

The PR model also has key predictions for loop pro- 
duction [l|, Q , and shows how the small-scale structure 
accounts for the production of loops at the small-scale 
cutoff in Nambu-Goto (NG) simulations. We were there- 
fore led to investigate loop distributions in AH simula- 
tions also, and to try to connect small scale structure and 
energy loss in our field theory simulations by looking for 
core (string width sized) loops. 

We found that the model for loop production in Ref. 
does not describe loop distributions in AH simulations, 
disagreeing with the observed power law. Furthermore, 
a lack of spatial correlation between core loops and long 
string, and a very low overall number density, roughly 
0.1 per horizon volume, point towards there being no 
significant energy loss through direct production of core 
loops from long string. 

Our results support the statement that significant en- 
ergy is lost directly from AH strings. This can be seen 
in snapshots from field theory simulations (Fig. [T|) with 
long string radiating energy directly into oscillations of 
the field as Higgs and gauge radiation. Our results also 
indicate that larger loops share this rapid energy-loss 
mechanism, which therefore needs to be added to the 
the standard model for loop distributions. We make a 
simple modification, that energy is radiated from loops 
at a constant rate. When coupled with an alternative 
model of loop production Q, which proposes fragmen- 
tation of horizon-sized loops, we are able to explain our 
results for the loop distribution. 

As an aside, a feature of the calculation for the loop dis- 
tribution function in Ref. is that the string correlators 
are assumed to be Gaussian, although it is argued that 
this is not crucial to the form of the result. The four point 
correlation in our simulations is explicitly non-Gaussian, 
but does scale. 

So when gravity can be neglected a general picture for 
the evolution of cosmic string networks is developing. A 
dominant energy loss mechanism works at the small-scale 
cutoff, which is loop production in NG simulations and 
classical radiation in the field theory. Energy is trans- 
ported from the horizon scale by small-scale structure. 
In both kinds of simulation loops are produced at the 
horizon scale, which then fragment with high probabil- 
ity. There is evidence of a population of non-intersecting 
loops remaining stable in NG [s^l, but in AH, where mas- 



sive radiation is available, loops radiate and shrink. 

We note that the small loops in NG simulations be- 
have like particles with mass ~ /i^i, which have similar 
average energy-momentum properties to classical mas- 
sive radiation. Thus the striking visual dissimilarity of 
NG and AH simulations is perhaps misleading: if one in- 
terprets the small loops as representing massive radiation 
and concentrates on the long string network, the differ- 
ences are greatly reduced. It is interesting to speculate 
that the small loops might really be massive states for 
fundamental cosmic strings. 

It remains a bit of a puzzle why AH strings have a 
highly suppressed production of small loops. The PR 
prediction, based on the NG equations of motion, is that 
small loop production should be dominated by the small- 
est physical scale, which for AH simulations is the string 
width, but few are observed. As the NG equations of 
motion break down at the string width scale, it is per- 
haps not surprising to see the break-down of a prediction 
based on those equations. 

In the traditional cosmic string scenario it is assumed 
that gravitational radiation back reaction plays a vital 
role in smoothing long strings, disguising the small-scale 
physics of the strings, and setting the scale for loop pro- 
duction. In this scenario most of the energy ends up as 
long-wavelength gravitational radiation, as opposed to 
ultra-high energy particles [13, [111 . 

It is by no means obvious that gravitational radiation 
reaction is strong enough to shut off the efficient and 
rapidly established energy loss mechanism afforded by 
massive radiation in a field theory. Our simulations show 
that it is not necessary to postulate gravitational radia- 
tion reaction to achieve scaling and a consistent picture 
of string network evolution. Nonetheless, it is still pos- 
sible that gravitational radiation could operate as tradi- 
tionally supposed. As a first step, we intend to test the 
impact of massless radiation on string network evolution 
with an investigation of small scale structure and loop 
distributions in field theories that have massless degrees 
of freedom, for instance global strings, and the Abelian 
Higgs model coupled to a dilaton. 
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